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In this work we consider the renormalization of the chiral two-pion exchange potential with 
explicit A-excitations for nucleon-nucleon scattering at next-to-leading (NLO) and next-to-next-to- 
leading order (N^LO). Due to the singular nature of the chiral potentials, correlations between 
different partial waves are generated. In particular we show that two-body scattering by a short dis- 
tance power like singular attractive interaction can be renormalized in all partial waves with a single 
counterterm, provided the singularities are identical. A parallel statement holds in the presence of 
tensor interactions when the eigenpotentials in the coupled channel problem also coincide. While 
this construction reduces the total number of counterterms to eleven in the case of nucleon-nucleon 
scattering with chiral two-pion exchange interactions with A degrees of freedom, the differences in 
the scattering phases as compared to the case with the uncorrelated partial wave renormalization 
become smaller as the angular momentum is increased in the elastic scattering region. 
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I. INTRODUCTION 



The basic and fundamental problem of nuclear physics 
is the determination of the nucleon-nucleon (NN) interac- 
tion [Ij . Field theoretical approaches to the nuclear force 
state that the NN potential can be expressed as a sum 
of increasingly heavy meson exchange contributions (for 
reviews see e.g. 0,3 ^-^id references therein). The result- 
ing high-quality potentials describe neutron-proton and 
proton-proton scattering data with a x^/d.o.f. < 1 0-0I- 
They all include at large distances the charge dependent 
one pion exchange (OPE) potential and typically need of 
the order of 40 parameters for parametrizing the shorter 
range components of the interaction. On the contrary, 
quantum chromodynamics (QCD), the underlying funda- 
mental theory of the strong interaction, requires only two 
parameters in the isospin symmetric limit: Aqcd and the 
average up and down quark masses. These fundamental 
QCD parameters can be traded for experimentally ac- 
cessible observables such as f-^, the pion weak decay con- 
stant, and m^r, the averaged pion mass. Obviously, the 
large number of parameters needed in phenomenologi- 
cal approaches arises within very specific schemes and 
functional forms. It is not clear whether this number of 
parameters can effectively be reduced by invoking rele- 
vant QCD features while maintaining the quality of the 
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description at the same time. Ultimately, lattice ab ini- 
tio calculations of the NN potential, for which incipient 
results already exist [8|, 0] , will eventually solve the prob- 
lem. 

In the present paper we deal with a situation where a 
reduction of parameters arises within the context of the 
renormalization of the chiral potentials deduced from the 
effective field theory (EFT) approach suggested by Wein- 
berg 1 1(1 two decades ago (for comprehensive reviews see 
Refs. [Ill, [l2l)- These chiral potentials turn out to be 
singular interactions which exhibit an inverse power law 
behaviour ~ l/r" at distances below the pion Comp- 
ton wavelength, niT^r <C 1. In case they are attractive 
the resulting amplitudes are sensitive to short distance 
physics and require renormalization. As will be shown, 
the renormalizability of singular attractive potentials can 
be translated into a mathematical short distance con- 
straint on the scattering amplitude for different partial 
waves. In the simplifying case of two-body scattering by a 
central attractive singular potential only one counterterm 
is needed in order to renormalize all the partial waves. 
In the more complex case of NN scattering, where spin, 
isospin dependence and tensor forces are present, this 
number can rise to eleven counterterms provided that 
certain conditions are met. We review below the EFT 
approach from the perspective of the number of param- 
eters in a way that our results can easily be displayed. 

The main appeal of the EFT idea in nuclear physics lies 
in the promise of a model independent approach where 
the long and short range contributions to observables 
can be disentangled with the aid of a sensible hierar- 
chy, eluding the ubiquitous problem of fine tunings. This 



feature is explicitly displayed through the introduction 
of counterterms in the effective Lagrangian, which en- 
code the underlying but unresolved short distance physics 
and are organized according to a power counting. The 
standard EFT formulation of the nuclear forces exploits 
the spontaneous breakdown of chiral symmetry, which re- 
quires derivative couplings for the pion. This implies that 
higher pion exchanges are power suppressed at momenta 
which are small compared with the chiral scale A^ ^ 
(47r/7r,Mjv) ^ IGeV. On the contrary, at high virtual 
momenta, pion exchanges become large, eventually re- 
quiring a suitable renormalization through the introduc- 
tion of counterterms. Furthermore, the non-perturbative 
nature of the NN interaction makes using a sensible re- 
summation of diagrams mandatory. A simple and effec- 
tive method already suggested by Weinberg [l^, [l^l and 
implemented for the first time by Ray, Ordonez and van 
Kolck [ij, [l^ proceeds in two steps: first, one deduces 
a chiral nuclear potential and afterwards solves the cor- 
responding Lippmann-Schwinger equation in momentum 
space, or equivalently, the Schrodinger equation in coor- 
dinate space. This scheme complies to the familiar and 
widely accepted concept of nuclear potential, which in 
the chiral case can be organized as an expansion in pow- 
ers of Q [H-il] 



V{r) = ^(0) W + V^^Hr) + V^'Hr) + V^^Hr) 



(1) 



where Q represents either the pion mass or the momen- 
tum of the nucleons (or additionally the nucleon-A split- 
ting in case we include the A-excitation as an explicit 
degree of freedom) . Within this scheme and making use 
of finite cut-offs at N'^LO [23,[23li the number of parame- 
ters becomes comparable with that of the phenomenolog- 
ical potentials even though the long distance behaviour 
is determined by chiral symmetry. 

A stringent constraint follows from the natural require- 
ment of short distance insensitivity: physics not explic- 
itly taken into account should be under control by fixing a 
sufficient amount of low energy parameters. Such a con- 
dition represents the basis of the renormalization process 
as understood in the present workQ. To fulfill this goal it 
is necessary to achieve approximate cut-off independence 
over a certain cut-off region. This immediately raises the 
question of what cut-off values can be regarded as nat- 
ural and how much the a priori arbitrary cut-off can be 



varied. The shortest de Broglie wavelength probed in 
elastic NN scattering below pion production threshold is 
Amin '^ 0.5fm, and thus we might expect stable results 
for similar short distance cut-offs. Otherwise the cut- 
off becomes an essential parameter of the theory. This 
particularly applies when the cut-off must be fine tuned 
to physical observables, a situation which actually takes 
place for specific power counting schemes. 

By naive power counting one expects chiral potentials 
to be singular for r <C l/m7r , 



V{r) 



1 






(2) 



where i' represents the order in the chiral expansion |f|, 
see Eq. ((IJ. While the resummations implied by solving 
the wave equation mix up the chiral power counting, they 
also enable finding non-perturbative new features when 
the short distance cut-off re becomes much smaller than 
any other long distance length scale. There arises the 
possibility of finding the adequate number of countert- 
erms which is compatible with the power counting of the 
long range potential and of obtaining a finite and unique 
limit for the corresponding scattering amplitude. In this 
regard several studies have found that the original Wein- 
berg power counting is inconsistent with renormalizabil- 
ity [2^,|3l,|3a]. This unexpected result has suggested sev- 
eral alternative approaches and heated debates question- 
ing the particular power counting, the renormalization 
process itself or the correctness of the non-perturbative 
resummation. We will not ponder on the pros and cons 
of any particular approach as this has been already done 
from several viewpoints in Refs. |29l . Isil |32| . ISa . |37| . At 
present it is unclear what aspects of the original EFT 
framework will ultimately prevail or be universally ac- 
cepted by the nuclear physics community (in this regard, 
see the related discussion in Ref. [38] ) . Rather than con- 
sidering the problem solved, we think that further work 
is still needed to settle the issue. 

In previous works [3^, S, I39l - t4l| we have shown that 
the necessary minimum number of counterterms renor- 
malizing a singular interaction can in fact be determined 
a priori from the behaviour of the configuration-space 
potential near the origin (for an earlier coordinate space 
treatment, see Ref. (42|). On a more numerical basis, 
similar observations have been made in momentum space 
either using finite cut-offs [29i, i41i| or subtractive meth- 
ods ■43-4(|. The main result is that in any uncoupled 



^ By renormalization we specifically mean the existence of well- 
defined scattering amplitudes when the cut-off is removed. This 
condition allows to identify all the short distance operators 
needed to remove the cut-off dependence, as done, for exam- 
ple, in Ref. 12911 for the OPE case. Once these counterterms are 
included in the computation approximate cut-off independence 
is assured and consequently there is no problem in keeping a fi- 
nite cut-off. For a different view on renormalization within an 
effective field theory context, see Lepage |30ll and the related 
discussions of Refs. |31|.|32||. 



^ The operator product expansion for six quark operators predict 
the functional form of the NN potential at short distance, which 
turns out to be a little weaker than 1/r^ and repulsive 133. |34|| . 
This shows that the NN potential computed on the lattice is reg- 
ular and hence might predict uniquely the NN scattering data as 
well as the deuteron properties and also that this short distance 
dependence is quite different from the chiral potentials. 

^ The inclusion of static degrees of freedom, such as the A isobar 
excitation in the small scale expansion, can change the expected 
power law behaviour of the potential. 



partial wave where an attractive singular potential is 
present, a counterterm is needed in order to renormal- 
ize this particular partial wave. An interesting corollary 
is that finite cut-off effects are less important the more 
singular the potential. Indeed, at order v in the chiral 



counting the potential behaves as 1/A 



i^+2„3+i/ 



generat- 



ing finite cut-off corrections (5a (fc) to the renormalized 
phase shifts 5ao{k) which scale as 

5Uk)-5K{k)^0{K-'''^--'^), (3) 

for large enough cut-offs [4Ci |41| , meaning in particular 
that cut-off independence is achieved in this case. On 
the contrary, if the potential is singular and repulsive 
the effect of counterterms becomes negligible for small 
enough cut-off radii. 

A particularly problematic consequence of the non- 
perturbative treatment of singular potentials is that the 
renormalization of two body scattering by an attractive 
singular central interaction requires an infinite number 
of counterterms, one per each partial wave. The un- 
limited proliferation of counterterms when renormaliz- 
ing singular interactions has been, among others, an ar- 
gument against removing the cut-off, as the resulting 
effective field theory will be unable to predict observ- 
ables [31|. This problem can be cured in perturbative 
power countings, like the one proposed by Kaplan, Sav- 
age and Wise 47, 48] , where the proliferation of countert- 
erms is naturally limited by the order of the approxima- 
tion. Unfortunately the singularity in the tensor forces 
makes the previo us p roposal poorly convergent in the 
'^Si — ^Di channel [43] (see however 50] for a renewed for- 
mulation) . In the modified Weinberg proposal of Nogga, 
Timmermans and van Kolck [20| , OPE is iterated in low 
angular momentum waves while treated perturbatively in 
sufficiently peripheral waves (usually I > 2). This choice 
naturally limits the number of necessary counterterms 
and, although it has has been criticized as arbitrary in 
Ref. 131,], it is sustained by the perturbative analysis of 
Ref. [37|. Higher order corrections are treated in pertur- 
bation theory and the corresponding (finite) number of 
counterterms is determined by imposing cut-off indepen- 
dence on the results, generating convergent amplitudes 
for the central waves with the NLO and N^LO chiral 
potentials [5l|. However, as we will show, there is a 
non-perturbative way of also obtaining a finite number 
of counterterms. 

In a recent paper [52| we have analyzed the role of A 
degrees of freedom for the central waves and the deuteron 
with the chiral A potentials of Refs. [ij, [Sj] with a rea- 
sonable phenomenological success. This particular po- 
tential furnishes simultaneously the theoretical require- 
ments of renormalizability □ and power counting. Actu- 
ally, convergence is achieved for reasonable cut-offs of the 



order of Tc ~ 0.5fm, that is, scales comparable with the 
shortest wavelength probed in NN elastic scattering be- 
low pion production threshold 0- Moreover, as discussed 
in Ref. [50], it may provide a good starting point for nu- 
clear matter calculations as it has a rather small _D-state 
probability, implying a sufficiently small wound integral 
which insures better convergence properties for the few 
body correlations and the nuclear many body problem. 

In the present work we analyze further the TPE po- 
tential with A-excitations [13, ISJ] , addressing the calcu- 
lation of non-central partial waves. We show how the 
number of counterterms can be made finite by imple- 
menting a renormalization prescription correlating an in- 
finite number of partial waves. Thus, this is a compelling 
example where, contrary to naive expectations, singular 
potentials may be consistently renormalized with a sin- 
gle common counterterm for all partial waves. The idea 
behind such a procedure is quite simple: if the potential 
has an inverse power law behaviour at short distances, 
V{r) ^ C/r^, with a coefficient C independent of energy 
and angular momentum, we expect all the reduced wave 
functions of the system to behave the same way at small 
enough radii, regardless of the energy or the angular mo- 
mentum, as the contribution from these two factors will 
become negligible in comparison with the strength of the 
potential. As a trivial consequence, all partial waves can 
be related to the zero energy s-wave. The issue is ana- 
lyzed in detail in Section |TT] both for regular and singular 
potentials. While this becomes a relevant observation for 
uncoupled channels, the tensor force requires a suitable 
generalization of the result for coupled channels, which 
is discussed in Sectioi i |IIII Surprisingly, the potentials 
computed in Refs. [13,1531 including TPE with A excita- 
tion do fulfill the necessary mathematical conditions to 
link partial waves with different angular momenta (see 
Section irvj) . Actually, we can estimate the finite cut-off 
error induced by these angular momentum correlations 
and which are exclusive of singular potentials. Other po- 
tentials do not automatically comply with these require- 
ments, so the question on the consistency of the partial 
waves correlations is not independent of the potential and 
indirectly on the power counting invoked to compute it. 
Of course, mathematical consistency does not necessarily 
mean phenomenological success, and we test our proposal 
against the widely accepted partial wave analysis (PWA) 
of the Nijmegen group in Section|Vl We see that actually 
there is no big difference between using the finite number 
of counterterms or renormalizing independently wave by 
wave, suggesting that improvements might be sought in 



* The divergence structure of this potential is identical to that 
of the chiral quark model in the Born-Oppenheimer approxima- 
tion \5^ , which being second order perturbation theory provides 
only attractive and singular potentials. 



^ This is a purely coordinate space argument where the cut-off 
in the potential has been removed . In momentum space the 
corresponding cut-off A ~ y'm^JJWjv ^ 350MeV implies also a 
regularization of the potential and an effective quenching of the 
dnNN coupling constant. This might be one of the reasons why 
momentum space calculations renormalizing the LS equation [4ll , 
[Sal require much larger cut-offs A ~ 1 — 4GeV than naively 
expected. 



the TPE chiral potentials as well. For completeness we 
address the problem of the familiar OPE in Appendix |X] 
. Finally, in Section IVII we summarize our main results 
and present our main conclusions and outlook for further 
work. 

As in previous works [33, lla, l39l - l41| we use extensively 
the coordinate space formulation which greatly simplifies 
the treatment and allows handy analytical calculations. 
In this approach, contact operators are treated implic- 
itly via boundary conditions in coordinate space. We 
do not follow an a priori power counting for the contact 
operators, but rather deduce the short range operator 
structure from the condition of cut-off independence and 
assuming that the long range potential is to be fully it- 
erated. The equivalence to momentum space renormal- 
ization has been discussed in detail for scattering states 
in Ref. [41| and the deuteron in Ref. |55| . 



the simplest possible contact interaction and it actually 
becomes equivalent to a short distance boundary condi- 
tion. We analyze below what can be obtained with such 
an interaction when the cut-off re is removed, both in the 
case of regular and singular interactions. 



A. Delta-Shell Potentials and Regular Interactions 

As mentioned, the two-body scattering problem can 
be described by the corresponding reduced Schrodinger 
equation, Eq. ^. For radii below the cut-off r^ there is 
no potential (due to the specific regularization employed 
for the finite range piece of the potential, see Eq. ([7])) 
and the solution for the wave function is simply 



Uk,i{i') — const X r 



for r < 



(9) 



II. CENTRAL DELTA-SHELL POTENTIALS 
AND THE PARTIAL WAVE EXPANSION 

As a preparation we will consider first the simplest two- 
body scattering problem described by a central potential 
V, which for the Z-wave reads 



'■k,l 



2nV{r) 



l{l + l) 



Uk.iir) ^k^Uk,i{r),{4^) 



where Uk^i is the reduced wave function, /i the reduced 
mass of the system, and k — y/TjlE is the center of 
mass momentum. The asymptotic long distance bound- 
ary condition is taken as 



Uk,i{r) 



kr 



Y + '^Kfc) 



(5) 



where Si{k) is the corresponding phase shift. We assume 
that V{r) can be decomposed as the sum of a finite range 
potential Vp (bounded by an exponential fall-off^ e"""") 
and a contact range interaction Vc 



V{r)^VF{r;r,) + Vcir;r,)^ 



(6) 



where we have added the auxiliary cut-off scale r^ which 
will be needed in order to regularize the contact range 
interaction. For convenience we have also regularized the 
finite range potential in the following way 



VFir;r,)^VFir)e{r-r,). 



(7) 



which means that the short range components of the fi- 
nite range potential are effectively absorbed in the con- 
tact potential Vc- For the contact potential we only con- 
sider for definiteness the case in which Vc is a delta-shell 
interaction 



Vcir;rc) 



47J-7-2 



■6{r - re) . 



(8) 



where Cq does not depend on energy, and no higher 
derivatives of the delta function are considered. This is 



where the regular solution has been chosen. 

The solution for radii above the cut-off depends (i) on 
the size of the cut-off with respect to the range of Vf 
and (ii) on whether Vf is a regular or singular interac- 
tion. For the present discussion, we will assume that the 
cut-off is much smaller than the range of Vf, which we 
will call gf, re <C a^, and that the finite range poten- 
tial is a regular one, limr^or^ Vpir) = 0. Under these 
conditions the reduced wave function can be written as 
a linear combination of a regular and irregular solution 
for r > Tf,, i.e. 



Uk,i{r) 



aiutfir) 



bi4:P(r) 



for r > re, (10) 



where the superscripts ^''''^^ and *^"''-' denote the regular 
and irregular solutions, respectively. 

For small radii, say re < r <C of, the behaviour of the 
regular and irregular wave functions is given by 



,(rog) 
^k,l 

(irr) 



(r) 
(r) 



.'+1 



(11) 

(12) 



where corrections depending on the presence of the po- 
tential Vf or the finite momentum k do not appear until 
higher relative powers of r are considered. 

The effect of the delta shell potential in the Schrodinger 
equation, Eq. ^, is to generate a discontinuity in the 
first derivative of the reduced wave function at r = Tj,. 
The previous statement can be summarized in the fol- 
lowing relation between the logarithmic derivatives of the 
wave functions for r < re and r > re 

2^ Coire) Mre) u^^f (re) + k{re) u'ifire) 



47rr? 



.(reg)/ 



.(irr)/ 



Mrc)u'Ci i^e) + hire) u'^}" ire) 

l + l 



(13) 



From this expression it can be checked that if we want 
the effect of the delta function to be nontrivial, we need 



Co{rc) to be a running coupling constant. In fact, for a 
constant value of Cq one finds in the r^, — > limit that 



ai{rc) 



^• 



,(rcg) 



(re) 



^2(+l 



(irr) 



(re) 



(14) 



meaning that the regular solution is effectively chosen as 
\bi\ *^ |«i|- Therefore, in order to avoid a trivial or irrel- 
evant contact interaction one needs that Co (re) evolves 
with Tc in a very specific way, a dependence that can be 
obtained by solving Eq. (fT3|) for a given bi/ai value. 

The running of Co (re) is so strongly determined by the 
scaling properties of the regular and irregular wave func- 
tions near the origin that, if Co(rc) is set to be non-trivial 
in a given partial wave, it will become trivial in all the 
other waves. This can be checked as follows. If we fit 
Co(rc) to reproduce bi/ai in the partial wave / = li and 
call this counterterm Cq (^c), its exact value can be ob- 
tained from solving Eq. ((T3)) for I = li. Using now the 

counterterm Cq ^ {re) for computing the linear combina- 
tion of solutions for / = hi^ h), we get the following 



bh 
ai2 



for I2 ^ h 



(15) 



which is just the same scaling as the corresponding one 
for a constant counterterm, Eq. (|14l) . Therefore we can 
take the simplification that a given counterterm only acts 
on a determined partial wave when the cut-off is removed, 
as is usually assumed. This means de facto a total inde- 
pendence of non-trivial counterterms for any other par- 
tial wave. Note that only a trivial counterterm produces 
a short distance interaction common to all partial waves. 



B. Delta Shell Potentials and Attractive Singular 
Interactions 



As we have seen, in order to have a non-trivial effect, 
the running of the counterterm depends on the scaling 
properties of the wave function near the origin, r — ^ 0. 
Given the fact that for regular potentials the scaling is 
different for each partial wave, the result is that only one 
partial wave will be affected by a given counterterm. On 
the contrary, as we will see, for attractive singular poten- 
tials the scaling does not depend on the angular momen- 
tum. Therefore, the scaling is independent of the partial 
wave chosen, and a given counterterm will affect all the 
partial waves simultaneously. This is our key observation 
which we will extend to tensor forces in Section IIIII and 
put forward below for the relevant case of TPE chiral NN 
interactions with A-excitations in Section IIVI 

Indeed, if we consider the behaviour of the reduced 
wave function for a power-law attractive singular poten- 
tial which for short enough distances behaves as 



2/iFF(r) 



m 



(16) 



where n > 2 and Rp is some given length scale which 
sets the range of the power-law behaviour of Vp- This 
new scale Rp may not coincide with the generic range 
ap of the potential as several lower energy scales may be 
present in the system |j. For a potential like the previous 
one and for distances Tc < r <^ Rp, the reduced wave 
function can be described by the WKB approximation 
since the de Broglie wavelength is slowly varying, 

Rp4- , ^ ^{r/RF}"^^-' « 1 , (17) 



yielding the short distance behaviour 



Uk.i{r) 





2 
n-2 


(Rf 


r Rp:^ r > r 


C 7 





1/2-1 



+Mk) 



(18) 



where Ai is some normalization constant and y/(A:) is a 
short distance phase which in principle depends on the 
angular momentum and the energy. For r < r^ the 
reduced wave function uj, / will show the expected r'+^ 
behaviour, see Eq. ©. 

Taking into account the behaviour of the wave func- 
tion around the cut-off, we can rewrite the equation that 
describes the running of Co for the case of singular inter- 
actions as 



2MCo(re) 
47rr? 



2 
R 

X cot 

l + l 



Rp 



n/2 



Rp 



•1/2-1 



Vi{k) 



(19) 



As can be immediately realized, for Tc — > the explicit 
/-dependent term stemming from the behaviour of the 
wave function for r < r^ can be dropped, leading to the 
following ultraviolet behaviour for Co 



2/iCo(re) 
47rr? 



-^ 



2 



R 



X cot 



1/2 



Rp 



i/2-l 



Vi{k) 



(20) 



which does not depend explicitly on angular momentum. 
In particular the previous equation means that for re — >■ 
we have the following identifications 

fi{ki) ^ (pi{k2) and ipi^{k) = Lpi^{k) . (21) 

In other words, the short distance phase is independent 
of angular momentum or energy. In the more general 



Actually in the chiral NN potential this is mostly the case where 
ap ~ 1/m-n- and Rp ~ ^/ f-K- 



case that we accept energy dependent counterterms Ck 
or higher derivatives of the delta, we can obtain an energy 
dependent semi-classical phase (pi{ki) 7^ (pi(k2), but the 
angular momentum independence will hold Q The only 
way in which one can break the condition ipi-^ (fc) — tpi^ (fc) 
is by accepting terms explicitly depending on the angular 
momentum in the contact range interaction. 

The previous results can be efSciently cast in the lan- 
guage of short distance boundary conditions as follows 



*fc,ii 



Uk,h 



^kM 



r— Tc-t-e 



'^k,U 



(22) 



r^r^+e 



that is, the logarithmic derivative at the cut-off radius 
of the reduced wave function does not depend on the 
angular momentum. This is the form in which we will 
effectively implement the condition ifii-^^k) — ipi^{k). 

It is important to notice that the angular momen- 
tum independence of the behaviour of the Z-wave reduced 
wave function near the origin is only realized when the 
cut-off is small enough, so the behaviour described in 
Eq. P^ is valid. In general Eq. P^ will be applicable 
for a radius such that the WKB approximation holds, 
and this radius will be smaller for higher partial waves. 
If we denote the previous radius as r^^^ , we will gen- 
erally have r^^^ > r^^^ for h < h- This means that 
for a given cut-off radius there will be a critical value of 
the angular momentum for which rj^^^ > Tc, and there- 
fore the condition of angular momentum independence 
should only be used for / < Ic, i.e. 



Uk,0 



-Tc+e 



Uk,l 



''k,l. 



-Tc+e 



Uk,lc-1 



(?3) 



-rc-t-e 



while for higher angular momenta the behaviour of the 
reduced wave function roughly corresponds to what is 
to be expected for a regular potential, and it may be 
possible to make a perturbative treatment as suggested 
inRefs. [sllll. 



III. THE INCLUSION OF THE TENSOR FORCE 

The previous analysis about the connection between 
the short range physics in different partial waves is only 
valid for the uncoupled channel case. If a tensor force 
is present, as in the nucleon-nucleon interaction, we will 



"^ The n-th derivative of the delta function potential generates a 
discontinuity on the (n+l)-th derivative of the wave function. As 
long as the behaviour of the wave function is given by Eq. III8I I , 
its (n + l)-th derivative will be angular momentum independent. 

* The ultraviolet or WKB behaviour described by Eq. I|18| l sets in 
once there is a deeply bound state in the system, or equivalently, 
once the zero energy wave function has reached its first zero 
(the finite energy wave function will contain anyway zeros due 
to the sin (fcr + 5i — Itt/2) behaviour at large distances) . This 



means that for radii bigger than r^ 
expansion is expected to converge. 



the perturbative 



need to account for the induced coupled channel structure 
happening in spin triplet channels in our analysis. For 
simplicity, we will assume a two-body system in which 
the finite range piece of the potential only contains a 
non-tensor (central) and a tensor piece 



Vpir^ = VNT{r) + Si2{f)VT{r) . 



(24) 



with 5*12 (r) = 3 (Ji • f(T2 • r — ffi -02 (i?i and 02 are the 
spin operators of particle 1 and 2), and where we do not 
specify any additional operator structure of the tensor 
and non-tensor pieces of the potential (like, for example, 
spin or isospin dependence) . The behaviour of the non- 
tensor and tensor piece at short enough distances is given 
by 



C 



NT 



VNTir) -^ 



Mr)^^, 



(25) 
(26) 



where we have assumed that they have the same power- 
law divergent behaviour at short distances. We have not 
yet determined whether the potentials are attractive or 
repulsive. 

The tensor force will couple spin triplet channels for 
which I = J ± 1, being the corresponding reduced 
Schrodinger equation for r > Vc 



- u" + UiiUj + Ui2Wj = k'^Uj , (27) 

~w'' + U12U.J + U22WJ = k^Wj . (28) 



where 



2i _ 2 id — 1) 

C/ii-2AilVT-2M7q^FT + ^^^^^, (29) 



2j + l 



Ui2=G 



Vj(J+i) 

2j + l 



2j; + 1 r^ 



(30) 



(31) 



The previous Schrodinger equation has four linearly inde- 
pendent solutions, of which only two of them are regular 
and therefore physically acceptable. The coupled chan- 
nel equations can be efficiently cast into the following 
compact notation 



u 



t2 



k'u 



where the wave function is now a matrix 



u 




(32) 



(33) 



with the *^°^ and ^^^ superscripts representing the two lin- 
early independent asymptotic (r — > 00) solutions of the 



system, and where the potential and the angular momen- 
tum operator are also 2x2 matrices: 



V = lVNT{r)+S{^VT{r), 

(j + l)(j+2) 



(34) 
(35) 



In the previous definition of the potential in matrix form, 
1 represents the identity and SJ2 the tensor operator, i.e. 



10\ 
1 ' 






1 



-2(j - 1) 6V7(J+T) 



2j + 1 K&VjU + 1) -2(J + 2) 



(36) 



(37) 



The tensor operator can be diagonalized with the follow- 
ing rotation 



rVo 



1 



%/7TT ^ 



VWTT \ -VI VITTJ ' 



leading to 






Lij j-j — rVj-L xv„ 



2 

-4 



j (j + 1) 2Vj(j + 1) 
2v/ja + l) i(j + 1) + 2 



where the d subscript indicates that the corresponding 
quantities are defined in the diagonal basis. In this ba- 
sis, the reduced Schrodinger equation is written in the 
following way 



- v" + 2m 


L^j, 2MCD(r,) 


= fc' Vj , 


where Vo = R 


., u," is the rotated wave fu 



5{7 



(41) 



A 



Rj V R represents the potential in the diagonal ba- 
sis, and the contact interaction has been explicitly in- 
troduced. At short enough distances, V^^d behaves in 
the following way 



1 I Cnt + 2 Ct 



V.D(.)->-(-"^0^^^ cJ-ACr)^ (42) 



A. Attractive-Attractive Case 

Situations where both eigenchannels are attractive are 
the easiest to handle. If a singular power-law potential 
is assumed, the behaviour of the reduced wave function 
for r < Tc can be safely ignored, yielding the following 
relation 



v/(?'c)vj(rc 



2pL 

Anr'^. 



Cu{rc), 



(43) 



meaning that there are three free parameters in this case 
(due to Cd being real and symmetric), in agreement with 
previous analysis of singular potentials in coupled chan- 
nels 36i, [SSf . The relation between the short range wave 
functions of channels with different total angular momen- 
tum is therefore 



T / T ' 



(44) 



(3o) where we have made use of Cd = C^. This relation is 
invariant with respect to the relative normalization of the 
two linearly independent solutions of each partial wave, 
which were previously denoted with the ^'^' and ^''' su- 
perscripts, and also on the set of linearly independent 

(39) solutions chosen, as can be easily checked. It should be 
noted too that the previous relation is reminiscent of the 
coupled-channel version of the two-potential formula of 
Ref. m. 

(40) A problem arises in relating the j — triplet state with 
other coupled channels, as the ^Pq wave is effectively an 
uncoupled state. 



Vi=o(r) 





V3p^{r) 



(45) 



Therefore, the previous representation of the short dis- 
tance potential, i.e. the coupled-channel delta-shell of 
Eq. (H5)) . cannot be a correct representation of the short 
range physics of the "^Pq channel. As a consequence the 
'^Pq channel cannot be unambiguously related with the 
other coupled channels. It is possible however to obtain 
convergent amplitudes by relating the ^Pq wave function 
with any of the "lower" eigen wave functions of the ref- 
erence wave function Vj , 



vip„{rc) _ (vj-)^i(re) 

V-^Poi^c) (Vj)2l(rc) 



vsp^jrc) _ i^jY22irc) 

V^Poi'^c) (Vj)22(?-c) ' 

(46) 



which means that depending on the signs and relative 
values of the van der Waals coefficients of the non-tensor 
and tensor piece, Cnt and Ct, the diagonalized poten- 
tial may be attractive-attractive, attractive-repulsive or 
repulsive-repulsive. Only two of these situations, namely 
the attractive-attractive and attractive-repulsive case, 
admit counterterms [36l . ISSj , and are therefore of interest 
from the point of view of renormalization. 



where 21 and 22 are the corresponding indices in the ro- 
tated wave function matrix W2- Both possibilities yield a 
renormalizable phase for the ^Pq wave, but not the same 
one: the phase shift depends on which of the previous 
two equations is used, indicating the presence of model 
dependence in any of these choices. Therefore, the only 
way to avoid model dependence is to treat the "^Pq wave 
as an independent wave in the attractive-attractive case. 
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B. Attractive- Repulsive Case 

If one of the eigenchannels is attractive and the other is 
repulsive, then we have that Eq. (P5| can only be applied 
in the attractive eigenchannel. In particular, the delta- 
shell coupling matrix Cd(?'c) takes the simplifying form 



CD(r-c 



'BC 



CA{rc)5AB5AC , 



(47) 



where the labels B,C — 1,2 are matrix indices, and A 
represents the index associated with the attractive solu- 
tion. That is, only one counterterm is needed in order 
to renormalize an attractive-repulsive coupled channel, 
as delta-shell contributions become trivial in the r-c — > 
limit, except if they happen in the AA subchannel. Com- 
bining the previous result with the boundary condition 
induced by the delta-shell potential 



j'irc) 



2^ 
47rr? 



Cu{rc)vjirc); 



(48) 



the following renormalization conditions are obtained 



(v,/(re)) 



AA 



i^n'irc)) 



AA 



i^3i'irc))AR 



(Vj2(rc))^^ 
i'^J2'irc))AR 



(vji(r-c))^fl^ (vi2(^c))^^ 

(Vjl'(rc) )iiA = (Vj2'(»-c) )iiA = ^ 

(vji'('^c) )iiii^ = i^n'irc) )rr^O., 



(49) 

(50) 

(51) 
(52) 



where R is the index of the repulsive eigenchannel. The 
first two equations relate the attractive eigenchannels 
of different partial waves, although they are redundant. 
From the form of the delta-shell coupling matrix Cd(?'c), 
together with Eq. (H51) . one obtains that 



i^/irc))AA _ i^Arc)) 
(vj(rc))^^ {vj{rc))AR 



AR 



(53) 



implying the equivalence of Eqs. (j49|) and ([50|) . The 
last two equations, (fSTj) and (f52l) . are just regularization 
conditions for the repulsive eigenchannels. Contrary to 
what happened in the attractive-attractive case, in the 
attractive- repulsive case the j = ^Pq wave (if attrac- 
tive) can be directly related with the other coupled triplet 
waves. 



IV. APPLICATION TO NEUTRON-PROTON 
SCATTERING 

The previous results can be applied to the case 
neutron-proton (up) scattering in nuclear effective field 
theory, where the resulting potentials are in many cases 
singular and attractive. In principle we can relate any 
two partial waves for which the potential diverges in the 
same way, i.e. the same power and the same coefficient. 
This means that correlations will emerge only between 
channels with the same spin and isospin values. For 



the spin triplet channels it is also necessary to consider 
whether the channels are coupled (Z = j ± 1) or uncou- 
pled {I — j). Therefore, we obtain a total of six sets of 
correlated waves, namely (i) singlet isovectors, (ii) singlet 
isoscalars, (iii) uncoupled triplet isovectors, (iv) uncou- 
pled triplet isoscalars, (v) coupled triplet isovectors, and 
(vi) coupled triplet isoscalars . We limit ourselves to the 
j < 5 partial waves. The "^Pq wave remains uncorrelated 
in our current scheme and it is not further considered in 
this work. 



A. The Chiral Neutron-Proton Potential 

The finite range piece of the nucleon-nucleon (NN) po- 
tential in chiral perturbation theory is expressed as an 
expansion in powers of Q 

VNN(r) = V^'^^r) + V^^Hf) + V^'\f) + 0{Q^) , (54) 

where Q represents either the pion mass, the nuclcon-A 
splitting or the momentum of the nucleons. We only con- 
sider here chiral potentials in which the A isobar has been 
explicitly included [13, [Sj] • The reasons for this decision 
are that (i) they have better convergence properties than 
their A-less counterparts, and (ii) they are attractive- 
attractive singular interactions in all coupled channels at 
orders Q^ (NLO) and Q^ (N^LO), leading to a simpler 
analysis in general. For the finite range piece of the inter- 
action we adopt the original Weinberg counting [lO|, in 
which I/Mat corrections are treated as higher order as it 
is done in Refs. [H, [60, |61|- At the orders considered in 
this work, the potential can be decomposed as a central, 
spin-spin and a tensor component which in coordinate 
space reads 

144N(r) - Vc{r) + TWc{r) + o {Vs{r) + TWs{r)) 

+ Si2{f){VT{r)+TWT{r)), (55) 

where spin-orbit and quadratic spin-orbit terms have 
been ignored as they do not appear up to higher or- 
ders. The use of previous counting rule for the l/M^ 
corrections is necessary if we plan to correlate at short 
distance the behaviour of the different waves, as it gen- 
erates a spin-orbit term which is less singular than the 
other components of the interaction. The operators t, a 
and 5*12 are given by 

T = fi-f2 = 2i(i+l)-3, 
Cr = (7l • (72 = 2s(s + 1) - 3 , 

Si2{r) = 3 CTi • rCT2 • r - CTi • (72 , (56) 

where Ti(2) and CTi(2) are the proton(neutron) isospin and 
spin operators; t and s represent the total isospin t = 0, 1 
and total spin s = 0, 1 of the np system. The precise form 
of the chiral A potential is taken from Ref. [53| . 

Note that in the singlet channel cases (s = 0) the ten- 
sor force operator does not contribute. For symmetry 
reasons (Fermi-Dirac statistics) we have (—)'+"+* — — l^ 
where I is the orbital angular momentum. This means in 



particular that even partial waves are isovectors (i = 1) 
and odd partial waves are isoscalars {t = 0). The NN 
potential reads for the singlet channels 

= Vc{r) + Wc{r) - 3Vs{r) - 3Ws{r) , 
Vip,{r)^Vip,{r)^Vi„,{r) 

= Vc{r) - 3Wc{r) - 3Vs{r) + 9Ws{r) , 

(57) 

that is, all the singlet channels can be described with two 
different potentials depending on whether we are in the 
isoscalar or isovector case (or equivalently, on whether 
even or odd partial waves are considered) . 

In the spin triplet channels (s = 1) we must distinguish 
between uncoupled (/ = j) and coupled waves (I = j±l). 
In the uncoupled waves, we can again distinguish between 
the potential in the isoscalar ('^-D2, ^^4) and isovector 
(^Pi, ^1)2) waves 

VsD.Ar) = V.GAr), 
VspAr)^VspAr)^V.HAr), 

(58) 

where the explicit expressions of the previous potentials 
in terms of the central, spin-spin and spin-tensor compo- 
nents is given by 

V.dA^) = Vc{r) - SWc{r) + Vsir) - SWs{r) 

+2VT{r)-6WT{r), (59) 

V.p^r) = Vc{r) + Wc{r) + Vs{r) + Ws{r) 
+2VTir) + 2WT(r) . 

(60) 

Equivalently, for the coupled waves we have 

Ri Vac, (r) Rj = R3 V3C3 (r) RJ = R5 V^c, (r) R 
R2V3c,(r)Rj=R4V3c,(r)Rl, 



T 

■5 I 



(61) 

for the isoscalar (^Ci, ^Cs, ^Cs) and isovector {^€2, ^€4) 
waves, and where the expressions are more cumbersome 
as they involve matrices. In the expressions above, Rj are 
the rotation matrices defined in Eq. ([55)) . The notation 
^Ci, 3(72, etc, is a short hand for ^Si-^Di, ^P2-^F2, etc, 
and the explicit form of the rotated potentials is given by 

Ri V3c,(r)RT = 1 (Vcir) - SWc{r) + Vsir) - SWs{r)) 

+ Si2,D iVAr) - SWrir)) , 
R2 V3c,(r) Rj = 1 (Vcir) + Wc{r) + Vsir) + Ws{r)) 



+ S{^^^{VT{r) + WTir)) 



(62) 



with 1 the 2x2 identity matrix, and Sj2 d the diagonal- 
ized tensor matrix represented by Eq. p9p . 



B. Van der Waals Behaviour of the Chiral 
Potentials 

At distances below the pion Compton wavelength, 
rnj^r <C 1, the chiral potentials exhibit at orders Q^ and 
Q^ a power-law behaviour of the type 



FM(r) 



C, 



M 



(63) 



'<{'^) 



with v — 2,Z and where the value of Cg depends on the 
particular component of the potential considered. These 
coefficients were computed in Ref. _52j based on the spec- 
tral function representation of the potentials of Krebs, 
Epelbaum and Meil3ner [5j|- It should be noted though 
that the exact behaviour of the potential at short dis- 
tances is inessential for the angular momentum correla- 
tions. What really matters is (i) that the potential is 
a singular attractive interaction and (ii) that it is much 
stronger than the centrifugal barrier at the chosen cut-off 
radius re- In the case of the order Q^ and Q^ chiral A- 
fuU potentials these conditions are fulfilled in all partial 
waves with j < 5 for cut-off radii as big as 1.0 fm. 



Correlated Renormalization of the Uncoupled 
Waves 



We describe the scattering states in the uncoupled 
waves by solving the following reduced Schrodinger equa- 
tion for r > rr 



'■k,l 



Mjv V^nn(7-) 



l{l + l) 



Uk,i{r) = k'^Uk,i{r) , 
(64) 



where u^^i is the reduced wave function, Vnn the corre- 
sponding chiral potential for the particular partial wave 
considered, k the center-of-mass momentum, / the an- 
gular momentum, and Mn is twice the neutron-proton 
reduced mass, i.e. Mn = 2MpMn/{Mp + M„). The re- 
duced wave function is asymptotically normalized to 

Uk.i (r) ^ fc' (cot Si ji (kr) - yi (fcr)) , (65) 

for r — > 00, with 5i the phase shift, and ji{x) = xji{x) 
and yi (x) = xyi {x) the reduced spherical Bessel func- 
tions. The normalization factor fc' is added in order to 
have a well-defined normalization of the wave function in 
the fc — ^ limit. At r = Tc the wave function can be 
determined by several means. One is by solving Eq. p^ 
for some value of the counterterm Co, which can be later 
fitted to reproduce some observable, like for example, the 
^5*0 (or ^Pi/^Pi/^D2) scattering length. A different way 
is to construct an asymptotic wave function (r — >■ 00) 
reproducing the desired scattering length. In the case of 
the ^^0 channel, this wave function is given by 



Wo,iSo(^) 



1 



r 
00 



(66) 
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with ao the ^5*0 scattering length, and then integrate the 
reduced Schrodinger equation, Eq. (j64l) , downwards from 
r — >■ oo to r = Tc. Then we use the different relations 
derived previously to obtain the logarithmic boundary 
condition at r = Vc for the different energies and partial 
waves considered. For the particular case of the ^5*0 chan- 
nel and its correlated channels ^D2 and ^G^, the relation 
takes the form 



-^0,150 



"O.iSo 



'■k^So 



^k^^So 



"fc,iSo 



"fc,iSo 



u 



k^D2 



"fc 1_D2 



^fc,lG4 



Uk.^Gi 



(67) 
(68) 



where the first equation relates the zero and finite energy 
states of the ^Sq wave, and the second one represents 
the partial wave correlations. For the other correlated 
channels, we have the correlation conditions 



^fc,iPi 



Wfc,iPi 



fe,3Pi 



''k.^Fa 



Uk.^Pi 



"^k^Fa 



k,^F3 



'■k.^H,, 



U 



k,^D2 



Uk,^D2 



Uk,^F3 



''k,^Gi 



^k,^Hr, 



'■k^Hs 



(69) 



"l^k^Hs 



Uk.^Gi 



(71) 



which are to be supplemented with the regularization 
conditions for the base waves 



(72) 
(73) 
(74) 



Mo 


IPl 


Uq 


IPl 


< 


^Pl 


Uq 


■^Pl 


K, 


^D2 



Uo,^D2 



"fc 


IPI 


Uk 

< 


IPI 


Uk 

< 


^Pl 

3P>2 


Uk 


^D2 



These boundary conditions can be used as initial integra- 
tion conditions for the corresponding Schrodinger equa- 
tion Eq. (IMl) . After integrating upwards from r = Tc to 
r — > oo, we match to the asymptotic behaviour of the 
wave functions, Eq. ([66]), in order to obtain the phase 
shifts. The equivalent value for the counterterm cou- 
pling Co (re) can be obtained from Eq. p^ . giving in the 
Tc -^ limit 



47rr^ 



^0,1 So 



(^c) 



Uo^Soirc) ' 



(75) 



plus the corresponding expressions for the other base 
waves. 



D. Correlated Renormalization of the Coupled 
Waves 



For the coupled channels we solve the coupled reduced 
Schrodinger equation in its matrix form 

L2' 



U/c. 



2/xVNN(r) 



Ufcj(r) = fc^Ufc.j(rj76) 



where we now use the notation of Sect. (|IIip in which u^ j , 
Vnn and L^ are matrices. The reduced wave function is 
asymptotically (r — > oo) normalized to 

Uk,j ^ (j, (fcr) M, (fc) - y, {kr)) F, (fc) , (77) 

where jj(fcr), yj{kr) and Fj{k) are diagonal matrices de- 
fined as 



yj-i{kr) 

Vj+i{kr) 

fcJ-i 
fc^+i 



yj{kr) 
F,(fc) 



(78) 
(79) 
(80) 



with ii{x) — X ji{x) and yi{x) = xyi{x) the reduced 
spherical Bessel functions. The normalization factor 
Fj(fc) is included in order to have a well-defined asymp- 
totic (r — > oo) wave function at fc = 0. The Mj(fc) 
matrix is the coupled-channel equivalent of cot 5i{k) and 
(70) is related to the S-matrix Sj(fc) by Mj(fc) = J(Sj(fc) + 



l)/(Sj(fc) — 1) with 1 the 2x2 identity matrix. 

For the chiral A-full potential of Ref. 153] all the cou- 
pled channels are attractive-attractive singular poten- 
tials at distances below the pion Compton wavelength. 
Thus three renormalization conditions or counterterms 
are needed in order to obtain well-defined results. The 
usual procedure is to fix the asymptotic (r — > oo) be- 
haviour of the wave function at fc = 0. That is, we fix 
the three scattering lengths of the coupled system. Then 
we integrate the Schrodinger equation, Eq. (j76p . down- 
wards from r — > oo to r = Tc. If we define Lfcj(r) as 

UAr)^<,A^)^k,r\r), (81) 

then, the finite energy solution is constructed from the 
following boundary condition at r = Tc 

Lfcj(r,)=Loj(r,). (82) 

The procedure for correlating the different partial 
waves considered is similar to the one employed in con- 
structing the finite energy solutions, the only difference 
being the rotation to the basis in which the tensor force 
is diagonal. For the two sets of correlated coupled chan- 
nels, those with j = 1,3,5 and those with j — 2,4, we 
have the boundary conditions 



Ri Lfc,i(?'c) Ri"^ = Ra Lfc,3(''c 



Ra 



1R4 



R5Lfc,5('' 



R5 , 

(83) 
(84) 



R2Lfe,2(?'c)R2 =R4Lfc,4(r 

from which the Mfc.j(fc) matrix (and the corresponding 
phase shifts) can be obtained. 
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FIG. 1. (Color online) (Upper panel) ^So, ^D2 and ^6*4 phase shifts computed from Eq. (|67p . using the ^Sq scattering length 
as an input parameter with a coordinate space cut-off Vc — 0.3 fm. (Lower panel) ^Pi, ^Fs and ^H^ phase shifts computed 
from Eq. H69p . using the ^Pi scattering length as an input parameter with a coordinate space cut-off r^ = 0.3 fm. The light 
blue band is generated by varying the cut-off radius within the range re = 0.6 — 0.8 fm. 



E. Cut-off Dependence of the Phase Shifts 

The cut-ofF dependence of the phase shifts in the corre- 
lated renormalization procedure can be easily estimated 
by making use of the renormalization group analysis of 
boundary condition regularization of Ref. [40| ■ For sim- 
plicity, we only consider in detail the uncoupled channel 
case. According to Ref. [40|, the cut-off dependence of 
the phase shift for an uncoupled channel is given by 



dSi{k;rc) 

dVr 



Mw VnnC^c) 



l{l + l) 



+ L'irc) + Lli(rc) ul,{r,) 



(85) 



where Si{k) is the phase shift, r^ the cut-off radius, and 
with Uk,i, Vnn, k, I and M^ as defined in Eq. (|64p. In the 
previous formula Lk^i{rc) is the logarithmic derivative of 
the Uk,i reduced wave function at the cut-off radius, i.e. 



Lk.iirc) 



Uk,i{rc) 



(86) 



If we are correlating the /q- and Z-waves, we have for the 
logarithmic derivatives at the cut-off radius Vc that 



LkoMi'^c) 



Lk.iirc); 



(87) 



where the partial wave Iq is taken to be the base wave, 
i.e. the wave for which we have fixed the value of the 
phase shift at fc = fco (or the scattering length if fco = 0). 
By taking into account that the reduced wave function 
Uko.ioir) obeys the following Schrodinger equation 



<o,/o('') 



^0(^0 



1) 



Uko,i„ir) 



Mat Vnn (r) 

it is trivial to check that the logarithmic boundary con- 
dition for fc == fcoi ^ = ^0 fulfills the differential equation 



Mtv VNN(r-c 



^0(^0 + 1) 



+ L'k„M,i^c) + Ll,oi^c)=0, (89) 

which is also the differential equation obeyed by Lk.i{rc). 
In particular, the previous means that the cut-off depen- 
dence of the phase shift simplifies to 

d6i{k;rc) il{l + l)-lo{k + l) 



dvc 



{e 



kl) 



ul.iirc) 



(90) 



For cut-off radii such that 2ra.,rrc <C 1, the behaviour 
of the wave functions will be determined by the van der 
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Waals piece of the interaction, i.e. m| i{rc) ^ r^, up to 
oscillations (see Eq. (UHl)) for the chiral NLO- and N^LO- 
A potentials of Ref. [53]. This implies that the cut-off 
dependence of the phase shifts can be approximated by 



5i„(fc,re)-(5i„(fc,0)(X-(fc2-fc2)^4^ 



(91) 



for / = Iq (that is, the base wave), and 

di{k,r,) - Si{k,0) (X [l{l + 1) - loilo + 1)] r2 , (92) 

for I ^ Iq. At the end of Section |Vl we illustrate these 
expectations for the chiral A-potentials of Ref. [53]. It 
should be noted that Eq. ([TO]) implies that the corre- 
lated renormalization procedure only generates converg- 
ing phase shifts if the potential Vnn is singular, as ex- 
pected from the discussion in Section |lTl The extension 
of the previous results to coupled channels is straightfor- 
ward and leads to the same conclusion and cut-off depen- 
dence as the uncoupled channel case. 

V. NUMERICAL RESULTS 

As we have shown we can relate the phase shifts in 
different partial waves using the short range relation de- 
scribed by Eqs. (I67ll7ip . We take in our numerical com- 
putations /^ = 92.4MeV, m^ = 138.03MeV, 2^„p = 
Mn = 2MpMn/{Mp + Mn) = 938.918MeV, gA = 1-29 in 
the OPE piece to account for the Goldberger-Treimann 
discrepancy and gA = 1.26 in the TPE piece of the po- 
tential. The discussion of the standard OPE potential 
corresponds to the attractive-repulsive case and is rele- 
gated to Appendix [Xj We discuss here the TPE chiral 



potential with A excitations as obtained from Ref. 153 ) 
(however with the spectral cut-off removed). For Ha, the 
chiral couplings ci , C3 and C4 and 6 = 63 -t- bs we take the 
values corresponding to "Fit 1" of Ref. [53] (see table I 
inside the previous reference). 

All the partial waves are renormalized at a cut-off ra- 
dius Tc — 0.3 fm. This cut-off is small enough by far: 
in most partial waves the phase shifts have already con- 
verged in the range Vc = 0.6—0.8 fm. Smaller cut-off radii 
are in principle possible, but require too much comput- 
ing time for the higher partial waves, while cut-off radii 
larger than 0.8 — 1.0 fm yield amplitudes which depend 
linearly on the cut-off for low partial waves. The ap- 
pearance of the first deeply bound state usually happens 
in the 0.5 — 1.0 fm region, the exact location depending 
on the particular partial wave considered. Even in the 
case of G- and H-waves there are usually between two 
and three deeply bound states at Tc = 0.3 fm. At these 
distances the wave functions are dominated by the van 
der Waals behaviour of the NLO - A and N^LO - A 
potentials, meaning that the correlated renormalization 
procedure is guaranteed to work. This does not imply 
however that the low energy phase shifts are dominated 
by the singular structure of the chiral potentials at dis- 
tances below the pion Compton wavelength. In fact, as 



will be commented in the following paragraphs, the re- 
sults for peripheral waves do not significantly differ from 
those computed in first order perturbation theory j53| . 
The explicit cut-off dependence of the phase shifts is dis- 
cussed in more detail for selected partial waves at the end 
of this section. 

In Fig. (HI), we show the results for the singlet waves. 
For the ^So-^D2-^Gi (^ Pi-^ F^-^ H5) correlation we have 
taken as input parameter the ^Sq {^Pi) scattering length 
from the Nijrnegen II potential ]5], which was com- 
puted in Ref. [62[ yielding the result 015^ ~ — 23.727fm 
(aipj = 2.797 fm'^). For the ^G^ phase in the isovector 
channels and the ^F^ wave in the isosinglet, the phase 



shifts do not differ much from those obtained in Ref. [53] 
in the Born approximation (in the previous reference only 
waves with I — 2,3,4 were considered). The ^H^ is also 
very similar to the phases obtained in Refs. [l6|, [13] for 
the NLO-A potential. These waves are also quite sim- 
ilar to those obtained in Ref. [361 ] by renormalizing the 
N'^LO potential for the A-less theory in a wave-by-wave 
basis, that is, by fixing the scattering lengths separately 
in each of the channels to their Nijmegen II values. In 
general, peripheral partial waves will not notice too much 
the inclusion of the two pion exchange interaction or the 
A excitation and will behave very similarly as in first 
order perturbation theory. In this regard, the partial 
wave correlation is useful mainly as a way to renormal- 
ize all the peripheral waves with a minimum number of 
counterterms, but not necessarily as a real correlation. 
The only wave in which it can be effectively noticed is 
in the ^D2 one, in which the ^5*0 — ^I?2 correlation pre- 
dicts a scattering length of ai/), = — 1.728 fm for the 
^D2 wave 13, to be compared with an optimal scattering 
length of ai£,2 = -1.686 fm^ for which the N^LO-A po- 
tential effectively reproduces the Nijmegen II results for 
^LAB < 150 MeV. The previous observation indicates 
the necessity of the specific ^D2 wave N'^LO counterterm 
in order to reproduce the results in this partial wave. The 
predicted value greatly differs from the one correspond- 
ing to the Nijmegen II or Reid93 potentials [gj, namely 
aiD2,Nijm2 = -1.389fm^ and mn, . Rein 9.3 = -1.377fm^, 
which were computed in Ref. [62| . This discrepancy is 
however common in most effective field theory computa- 
tions in which the scattering length is fixed, sec for ex- 
ample Ref. [36], or the related comments in Refs. [iil - lia ). 
where a subtractive regularization approach is employed. 
This inconsistency between the Nijmegen low energy pa- 



^ The quoted scattering lengths have been computed for the fixed 
cut-off radius r^ = 0.3 fm for the N^LO-A case and are accurate 
within the numerical error. The systematic uncertainty of taking 
a cut-off between 0.3 and 0.6 fm typically only influences the 
last digit (in the particular case of the ^D2 scattering length we 
have a 0.005 fm"^ change in the previous cut-off window). Note 
that a factor of two in the cut-off corresponds to doubling the 
momentum space cut-off. Actually, using the "equivalence" A = 
7r/2rc |4ll. I55II . we are testing the A = 0.5 — 1 GeV region which 
seems reasonable. 
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rameters and the chiral A potentials is explained by the 
fact that the Nijmegen phenomenological potentials do 
not contain either two pion exchange contributions or A 
excitations, and was briefly commented in the previous 
references. 

In Fig. ^, we show the results for the uncou- 
pled triplet waves. We have taken as input parame- 
ter for the ^ Pi (^ Dy,) correlation the Nijmegen II scat- 
tering length [g^l, namely asp^ ~ 1.529 fm^ {o,'^D2 — 
— 7.405 fm ). Taking these values does not yield the bet- 
ter possible results for the ^Pi ('^-Di) wave, but generates 
renormalized results for the ^F^ — ^H^ {^G4J waves. As 
it happened in the singlet case, the phase shifts for the 
higher partial waves are very similar to the values ob- 
tained in first order perturbation theory either in the 
A-less [H and A-full [H [H cases. 

In Figs. ([3]) and ^ we show the results for the cou- 
pled triplet waves. For the '^5'i — ^Di correlation, we 
have taken as input parameters the deuteron binding en- 
ergy Bd = 2.224575 MeV and D/S asymptotic ratio r] — 
0.0256, and the '^Si scattering length aag.^ = 5.419 fm. 
The scattering solutions are then constructed by orthogo- 
nality with respect to the deuteron wave function and the 
"^Si scattering state. The procedure is described in detail 
in Ref . [33] , and was already used in Ref . [52] to construct 
the scattering solutions in the '^S'l — "^Di channel for the 
A potentials of Ref. p^. For re — 0.3 fm, we obtain 
the values aE^ — 1.953 fm and as^i^ — 5.034 fm for the 
scattering lengths. The previous low energy information 
yields much better results than the use of the Nijmegen 
scattering length for this channels (gei — 1.647 fm'^ and 
as El = 6.505 fm ), which induce a spurious resonance at 
fccm — 100 MeV when the N^LO-A potentials are em- 
ployed. This behaviour can also happen when using the 
standard chiral potentials without explicit A degrees of 
freedom, as was discussed in Refs. [44 - 46] . The corre- 
sponding phase shift for the ^D^ is slightly better than 
the one obtained in Ref. [5j], while the '^G'3 phase and 
the 63 mixing parameter are quite similar to the ones 
obtained in the previous reference. The results for the 
•^Ga — ^/a coupled channel are good in general with the 
exception of the '^G'5 phase in which only the threshold 
behaviour is correctly reproduced. 

In the case of the '^^2 — ^-^2 waves, we renormalize 
these waves by fixing the scattering lengths to the val- 
ues asp^ = — 0.320 fm , a^^ = 1.936 fm and asp^ = 
— 1.289 fm , which provide an acceptable description of 
the phase shifts for this channel, see Fig. (g]). As hap- 
pened in the ^Si — "^Di channel, the Nijmegen II values 
for the scattering lengths [62| do not yield good results 
with the NLO — A and N^LO — A potentials. However, 
in the A-less theory the Nijmegen II scattering lengths 
generated good results at N^LO [3a], meaning that the 
discrepancy is due to the long range physics introduced 
by the A excitations. The '^i^2 phases show a strong rela- 
tive cut-off dependence in the range re = 0.6 — 0.8 fm, al- 
though this is partly due to the small value of this phase. 
As in the previous cases, the ^F^ — ^H^ waves are very 



similar to those of Ref. '53] . 

Finally, in Fig. ^ we consider the cut-off dependence 
of the phase shifts in the form of a logarithmic (or Lep- 
age's [20]) plot. For simplicity, we only consider two sets 
of correlated partial waves in detail, ^5*0-^-02-^^4 and 
'^Pi-^F^-^H^, a singlet and a triplet. For the other par- 
tial waves the cut-off dependence follows a similar pat- 
tern. According to Sect. IIVEI for small enough cut-offs 
(2TO,n-rc ^ 1) the convergence of the phase shift in the 
lower partial wave of the correlation (i.e. ^Sq/^Pi in this 
case) is given by 

log \5A{k, re) - SAik, 0)] ~ 4 logr^ + Ca + fA{rc) , 

(93) 

where A = ^S()(^Pi), 6A{k,rc) is the phase shift com- 
puted at the cut-off radius re, SAik, 0) the phase shift in 
the Tc — >■ limit, Ca a constant, and /a (a;) a small os- 
cillatory contribution which takes into account the sine 
factor of the reduced wave function at short distances, 
see Eq. ^TE\\ . As can be seen in Fig. ([SJ, this behaviour is 
indeed fulfilled up to re ^ 0.8 — 0.9 fm. In particular, the 
numerical factors multiplying the logarithms in the fits 
of Fig. ([5]) are very close to 4, indicating that the van der 
Waals contribution to the chiral potential dominates the 
behaviour of the wave functions at short distances. For 
the higher partial waves in the correlation, the expected 
scaling is 



logl<5B(fc,r,)-fe(fc,0)l 



2 \ogre + Cb + fsire) , 

(94) 



with B = ^D2/^Gi{^Fz/^H^). In these waves the van 
der Waals dominance is apparent for cut-off radii below 
Tc ~ 0.5 — 0.8 fm, with the lower bound corresponding 
to the most peripheral partial waves. It should be noted 
however that the appearance of van der Waals scaling in 
the renormalization group (RG) flow of the phase shifts 
for the higher partial waves does not imply that the phase 
shifts themselves are dominated by the 1/r^ piece of the 
interaction. The region in which the RG flow is driven 
by the chiral van der Waals force only amounts for a tiny 
contribution to the total phase shifts of the peripheral 
waves, as can be deduced from the large negative values 
of \og\5{k,rc) — 5{k,Q)\ in the case of the ^6*4 and ^H<^ 
waves, see Fig. ([5]). This feature fully agrees with the 
expectations of the renormalization approach. 



VI. CONCLUSIONS 

In the present paper we have considered the relation 
between the renormalization of attractive singular poten- 
tials and the partial wave expansion. Given that attrac- 
tive singular interactions can be renormalized by includ- 
ing one counterterm per partial wave, each counterterm 
stabilizes the cut-off dependence in each one of the chan- 
nels separately. While this is a sufficient condition for 
renormalizability it is actually not necessary. We have 
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FIG. 2. (Color online) (Upper panel) ^Pi, '^F^ and ^H^^ phase shifts computed from Eq. H70p . using the ^Pi scattering length 
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Eq. (|7ip . using the '^D2 scattering length as an input parameter with a coordinate space cut-off r^ = 0.3 fm. The light blue 
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shown that if the finite range (attractive singular) inter- 
action is central, then it can be renormahzed by means of 
a single delta-shell central potential in coordinate space, 
in contrast with the previous situation in which the pre- 
dictive power is lost as there are an infinite number of 
partial waves. Of course, this result depends on the as- 
sumption that the unknown short range potential which 
is represented by a single delta-shell counterterm is cen- 
tral. Phenomenological potentials do depend on the or- 
bital angular momentum at short distances 0-01 ■ For 
this more general situation in which nothing can be as- 
sumed about the short range interaction, the usual result 
of one counterterm per channel will be recovered. 

Our analysis has been carried out in coordinate space, 
which on the other hand has been proven to be equivalent 
to momentum space calculations [4l|, [Sa] . The particu- 
larly interesting issue of extending the correlated renor- 
malization method to momentum space, not addressed 
in the present work, is left for future research. A possi- 
ble clue might be provided by the observation that the 
high momentum behaviour of the chiral potentials ought 
to reflect the partial wave independence observed and 
exploited in the present paper at short distances, sug- 
gestinga common subtraction perhaps along the lines of 
Refs. [iS-Iil. 

We have extended the previous result to the case of a 



finite range potential containing a tensor piece, which is 
of great interest for the renormalization of nuclear forces 
in the effective field theory approach. In that case, the 
number of counterterms depends on the sign of the eigen- 
values of the coupled channel potential. The application 
to the chiral NN potentials with A excitations is possi- 
ble and straightforward, and only requires to take into 
account the additional spin and isospin structure of the 
NN interaction. We stress that this is based on taking 
a counterterm structure based on the longer range OPE 
and TPE components of the interaction. For the order 
Q^ and Q^ A-potentials of Ref . [53.] a total of eleven coun- 
terterm is found to be needed to completely renormalize 
the interaction in all channels. This is only two more 
counterterms than what Weinberg's dimensional power 
counting dictates for the contact range interaction at the 
considered orders. 
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Appendix A: Partial Wave Correlations with the 
One Pion Exchange Potential 



In this appendix we review the correlated renormal- 
ization for the one pion exchange potential case, which 
corresponds to the leading order piece of the chiral po- 
tential. The OPE potential can be decomposed as 



VoPEir^=aTWs{r) + S,2ir)TWTir), (Al) 
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FIG. 4. (Color online) ^Pa 



- '^F2 and ^F4 



' ^^4 phase shifts. The ^Pa 



_F2 wave is constructed to reproduce the Nijmegen II 



scattering lengths for this coupled channel, while the Fa — Ha wave is obtained from the partial wave correlation described 
in Eq. (I84p without introducing new counterterms. The values of the coordinate space cut-offs are the same as in Figs. ^ and 



where the operators a, r and 6*12 were defined in Eq. (j56p 
and Ws (r) and Wt (r) are given by 



Ws{r) 
Wrir) 



487r/2 
487r/2 



1 



mi^r 



{m-^ry 



(A2) 
(A3) 



As can be seen, the only singular component of the OPE 
potential is the tensor piece. Therefore partial wave 
correlations only arise between attractive triplet partial 
waves. Specifically, there are three sets of correlated 
waves: (i) the ^Ci — "^Ca — ^C^ case, which happens 
between coupled waves, (u) the '^Pq — '^€2 — ^Ca case. 



in which there is one uncoupled wave (the '^Pq) and the 
rest are coupled, and (iii) ■^£'2 ^ ^^4 in which all waves 
are uncoupled triplets. All the coupled waves are of the 
attractive-repulsive type, and in total only three coun- 
terterms are needed in order to obtain finite scattering 
amplitudes for the OPE potential. Nonetheless it should 
be noted that in usual EFT computations a fourth coun- 
terterm will be added to renormalize the ^5*0 wave. In 
any case, we will only consider those waves which can be 
related. 

Contrary to the NLO - A and N^LO - A cases, the 
cut-off radius must be quite small in order for the partial 
wave correlations to converge (specially between the ■^Pq 
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and '^^2 — ^-^2 waves). In particular we take re — 0.15 fm. 
The ^Ci, ^Ca and ^Cs correlation is shown in Fig. ([6|). 
For this case, the ^Si — ^Di wave function is renormalized 
by reproducing the triplet ^^i scattering length, 035^ ~ 
5.419 fm, a procedure which was described in detail in 
Ref. [33] ■ The other partial waves are generated by the 
renormalization conditions given in Eqs. (I49ll52p . As can 
be seen, the description of the Ei and ^D^ wave is not 
especially good; these waves improve noticeably with the 



inclusion of two-pion exchange and the A. The remaining 
j — 3 and j = 5 phases do not differ two much from their 
NLO-A and N^LO-A counterparts, as expected from the 
fact that peripheral waves are OPE dominated. 

In Fig. d?]) we show the resulting ^Pg, '^^2 — ^^2 and 
•^^4 — ^H4 phase shifts, which have been obtained by 
using the ^Pq wave as the base wave, and where the "^Pq 
scattering length has been taken to be aip^ = —2.71 fm . 
We can see that the ■^Fo, E2 and '^^4 waves are not well 
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The coordinate space cut-off is taken to be Tc = 0.15 fm. 



reproduced with OPE alone and need the inclusion of the 
higher orders of the potential. 

Finally, in Fig. (j8]) the '^1)2 and '^G4 phase shifts are 
shown. The '^^2 phase has been renormalized to re- 
produce the Nijmegen II value of the scattering length, 



a^jj^ ^ -7.405 fm^ The OPE results for the ^D2 are 
worse than those of NLO-A and N^LO-A at moderate 
energies of the order of ii^LAB > 150 MeV. The '^G4 phase 
is nicely reproduced with OPE alone. 
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